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Abstract. We report on the recent developments of a new effective field theory for nuclear matter 
[1,2,3]. We present first the nuclear matter chiral power counting that takes into account both short- 
and long-range inter-nucleon interactions. It also identifies non-perturbative strings of diagrams, 
related to the iteration of nucleon-nucleon interactions, which have to be re-summed. The methods 
of unitary chiral perturbation theory has been shown to be a useful tool in order to perform those 
resummations. Results up to next-to-leading order for the ground state energy per particle of nuclear 
matter, the in-medium chiral quark condensate and pion self-energy are discussed. 
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INTRODUCTION 

An interesting achievement in nuclear physics would be the calculation of atomic nuclei 
and nuclear matter properties from microscopic inter-nucleon forces in a systematic and 
controlled way. This is a non-perturbative problem involving the strong interactions. In 
the last decades, effective field theory (EFT) has been proven to be an indispensable 
tool to accomplish such an ambitious goal. In this work we employ chiral effective field 
theory (^EFT) [4, 5] to nuclear systems [6, 7], with nucleons and pions as the pertinent 
degrees of freedom. For the lightest nuclear systems with two, three and four nucleons 
#EFT has been successfully applied [8]. For heavier nuclei one common procedure is 
to employ the chiral nucleon-nucleon (NN) potential derived in ^EFT combined with 
standard many -body methods. In [1] we have derived a chiral power counting in nuclear 
matter that takes into account local multi-nucleon interactions simultaneously to the 
pion-nucleon interactions. Many present applications to nuclei and nuclear matter only 
consider meson-baryon chiral Lagrangians without constraints from free NN scattering. 
Others in fact include short range interactions, but loose contact to the vacuum by fitting 
free parameters to nuclear matter properties. In [2] the techniques of unitary chiral 
perturbation theory (U^PT) [9, 10] have been applied to perform the resummation of 
the non-perturbative nature of the NN-interaction to vacuum phase shifts as well as to 
in-medium NN-scattering up to next-to-leading order (NLO). Our novel power counting 
has been applied in [1, 2, 3] to a variety of problems in nuclear matter. 



IN-MEDIUM CHIRAL POWER COUNTING 



We consider the ground state of nuclear matter which, under the action of any time- 
dependent operator at asymptotic times, behaves as Fermi seas of nucleons [11]. The 
Fermi seas are filled states up to a limiting Fermi momentum ^ = (37r 2 p^) 1 / 3 , with 
pi the density of the corresponding Fermi sea and A refers to the third component 
of the nucleon isospin. The nucleon propagator in the medium contains a particle part 
propagating above the Fermi surface and a hole part propagating below [12]. 

(1) 

In [11] there has been established the concept of an "in-medium generalized vertex" 
(IGV). Such type of vertices results because one can connect several bilinear vacuum 
vertices through the exchange of baryon propagators with the flow of a baryon number 
of one through a closed loop. At least one Fermi sea insertion is needed, otherwise we 
would have a vacuum closed nucleon loop that in a low energy effective field theory 
is buried in the higher order counterterms. In order to treat chiral Lagrangians with an 
arbitrary number of baryon fields (bilinear, quartic, etc.) we consider bilinear vertices 
like in [11, 13], but allow the additional exchange of any type of heavy meson fields. 
Those should be considered merely as auxiliary fields that allow one to find a tractable 
representation of the multi-nucleon interactions that result when the masses of the heavy 
mesons tend to infinity. Consequently a heavy meson propagator is counted as <ff(p Q ) 
due to their large masses. In addition, we employ the non-standard counting case from 
the start [13] and any nucleon propagator is considered as <f?(p~ 2 ). In this way, no 
diagram whose chiral order is actually lower than expected by following the standard 
counting rules is lost. We regard m K ~ ^ ~ &(p), which are considered much smaller 
than a hadronic scale A^. of several hundred MeV that results by integrating out all 
other particle types, including nucleons with larger three-momentum, heavy mesons and 
nucleon isobars [7]. The final formula for the chiral order p v of a given diagram is [1] 
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where E n ,V n , Vg and V p is the number of external pion lines, mesonic vertices, bilinear 
baryon vertices and IGVs, respectively. The number of pion lines attached to the i th 
mesonic vertex is indicated by n ; and ii corresponds to the chiral order of the latter, dj 
is the chiral order of the i th bilinear baryon vertex, v ; is the number of all mesonic lines 
attached to it and W[ the number of only the heavy lines. It is worth stressing that v is 
bound from below and eq.(2) implies conditions for augmenting the number of lines in 
a diagram without increasing the chiral power (see [1] for details). Because of the last 
term in eq.(2) adding a new IGV to a connected diagram increases the counting at least 
by one unit. The number v given in eq.(2) represents a lower limit for the actual chiral 
power fi of a diagram, because in some instances (like standard ^PT subdiagrams) the 
nucleon propagator may follow the standard counting. 



For evaluations we employ the chiral effective Langrangian 
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Xeff = £ ^ + £ + £ > (3) 
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where, from left to right one has the mesonic ^PT Lagrangian, the Lagrangian of heavy 
baryon %PT (HB^PT) [14] and the NN Lagrangian [6, 7]. The scattering amplitudes 
that follow these Lagrangians can be found in [2]. 



NON-PERTURBATIVE METHODS 



Figure 1. Resummation of the two-nucleon reducible diagrams. A wiggly line represents the 
interaction kernel Nji of fig.2. 

In [6, 7] it was argued that due to the large nucleon mass the two-nucleon reducible 
diagrams should be resummed, as it is schematically depicted in fig. 1 . This is also a 
consequence of eq.(2) when considering lowest-order NN-interactions. The wiggly line 
is understood as a local NN interaction plus the exchange of a pion as depicted on the 
left hand side of fig.2. 
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Figure 2. The exchange of a wiggly line between two nucleons indicates the sum of the local 
and the one -pion exchange. A one-wiggle exchange corresponds to the leading order (LO), Lyy, 
and the two-wiggle exchange box diagram to NLO, Ljj. 

We follow the techniques of U^PT [2, 9, 10] that performs the resummation partial 
wave by partial wave. The relation for the scattering amplitude, Tji, in U^PT is 



Tjj{W) = [l+Nji(£,e,S) -g] 1 -Nji{W) . 



(4) 



It results by performing a once subtracted dispersion relation of the inverse of a partial 
wave amplitude taking an integration contour, denoted by C/UC// in fig. 3, closed by 
a circle with infinite radius that engulfs both the left- and the right-hand cut. The 




Figure 3. Right- and left-hand cuts of T n (£,£,S). 



quantum numbers are total isospin /, total spin S, total angular momentum / and the 
out- and in-going orbital angular momenta i, £, in order. The function Njj(A) is called 
the interaction kernel. The unitarity loop function g(A) shares the right-hand cut of the 
scattering amplitude and its discontinuity along this cut is — im\k\/(2n). In this way, 
the effects of the nucleon mass associated with the two-nucleon reducible diagrams are 
taken into account. A once- subtracted dispersion relation can be written down given the 
degree of divergence for k 2 — > °°. The integration contour is shown in fig. 3 by Q. It 
results 
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S (A) = m - j g dk 2___ =,„-,— . (5, 

One subtraction has been taken at D < so that the integral is convergent, yielding 
the subtraction constant go. Note that g(A) is real along the left-hand cut. According to 
eq.(4) we can relate the imaginary part of the interaction kernel to the imaginary part of 
the scattering amplitude along the left-hand cut 
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We employ this result to write down a once- subtracted dispersion relation for Njj. The 
integration contour is shown in fig.3 as Qj and we obtain 

NjI (A) = Njj(D) + — J__ dk {k 2_ A _ ie){k2 _ D) • (7) 

This is a non-linear integral equation in Njj (A). However, one can obtain an approximate 
algebraic solution by choosing appropriately the subtraction constant go in eq.(5) so 
that g(k 2 ) vanishes at one point in the low energy part of the left-hand cut. Making 
g(—m 2 ) = so that \k\ = im K one gets go ~ — = —0.54 m 2 , the natural size for go- 
In this way, with g(k 2 ) suppressed in the low energy part of the left-hand cut, one can 
obtain a solution of eq.(7) as an expansion in powers of g. For practical terms this is done 
as follows. In a plain chiral calculation of NN partial waves the cuts are not resummed. 
But we can match the chiral calculation with an expansion of the geometric series in 
eq.(4) 

T n (£, I S)=N JI (£, £, S) - N n (£, £, S)-g-N JI (£,£,S) + ... (8) 

up to the same number of two-nucleon reducible loops. In addition, one has the chiral 
expansion of the interaction kernel 
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Nji = £ = NP + G{p 2 ) , (9) 
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where the chiral order is indicated by the superscript. In order to extract Njj, the <ff(p n ) 
chiral calculation of a NN partial wave comprising at most n two-nucleon reducible 
loops is matched with eq.(8) 
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On the right-hand-side of eq.(10) the chiral calculation is given, which is diagrammati- 
cally shown in fig. 2. From eq.(10) we obtain N$ = and 

<=4! , +<^< ) . (id 

which corresponds to the difference between a full calculation of a two-nucleon re- 
ducible loop and the result obtained by factorizing the vertices on-shell, so that the 
right-hand cut is removed. This is the reason why a two-nucleon reducible loop cal- 
culated in ^PT is counted as 6{p) in order to evaluate Njj. This fact is incorporated in 
the interaction kernel Njj, which can be improved order by order. 

The evaluation of the NN scattering amplitudes in the nuclear medium can be obtained 
from the results in vacuum, since the only modification without increasing the chiral 
order corresponds to use the full in-medium nucleon propagator G, eq.(l). This is 
directly accomplished by replacing 

g = L\Qj — > Lio = Li ,/ + 2Lio, m +L 10;d =Lio,p P +L 10) / ! / ! (12) 

in eq.(4), for any order of the NN scattering in the nuclear medium. The subscripts /, m 
and d indicate zero, one or two Fermi sea insertions from the in-medium nucleon propa- 
gators. The same process as discussed previously for the vacuum case is followed to fix 
Nji, only that in all loop integrals the free nucleon propagators have to be exchanged by 
their full in-medium equivalents. 



NUCLEAR MATTER ENERGY PER PARTICLE 
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Figure 4. Contributions to the nuclear matter energy density £ = (£l\J%ff\£l) up to 0{p 6 ). 
Pions are depicted by dashed lines. 

The diagrams required to study the problem of the nuclear matter ground state energy 
per particle E/A = S'/p, by applying eq.(2) up to NLO, are shown in fig. 4. The kinetic 
energy closed by a Fermi sea insertion is depicted by diagram 1 and corresponds to the 
energy of a free Fermi gas. All other diagrams are obtained by considering nucleon self- 
energy contributions which are closed with the full in-medium propagator G, eq.(l). 
depends on go both implicitly, due to the dependence of the partial wave amplitudes Tjj 
on go, and explicitly in a linear manner. For the latter dependence the symbol go was 



introduced in [2] (to which we refer for a more detailed discussion on go and go). For 
pure neutron matter both go and go are obtained with the value go = go — —0.6 m|, 
very close to their expected natural size go — —0.54 m|. The case of symmetric nuclear 
matter requires some fine-tuning of go, with a fi na l value of go — — m|, while go is kept 
on its natural value. Results in perfect agreement with experiments were obtained for 




Figure 5. Energy per particle E/A = $ /p vs. the density p for neutron (left panel) and sym- 
metric nuclear matter (right panel). 

the saturation density, energy per particle and compression modulus [2]. The equations 
of state for neutron (left panel) and symmetric nuclear matter (right panel) are shown in 
fig.5. The magenta lines are from the many-body calculations of [15] employing to so- 
called realistic NN potentials. The agreement between these references and our results 
is remarkable. 

IN-MEDIUM CHIRAL QUARK CONDENSATE 
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Figure 6. Contributions S, to the in-medium chiral quark condensate (£l\qfqj\£l) = 
~8/8sjj(x)^(v,a } s,p)\ VM . s _ p= oup to 0(p 6 ). The external scalar source in indicated by the wavy 
line. 

The diagrams contributing to the problem of the in-medium chiral quark condensate, 
by applying eq.(2) up to NLO, are shown in fig. 6. Since there are no explicit quark 



degrees of freedom incorporated in the theory, the scalar chiral quark condensate is 
emulated by a scalar source [5]. Diagram 1 can be identified with the sigma-term of 
pion-nucleon scattering plus a small isospin breaking contribution. Diagram 2 vanishes, 
while diagram 3 is suppressed by one order. In diagram 6 the scalar source couples to a 
one-pion exchange of the NN interactions. For the diagrams 4 and 5 the scalar source is 
attached to the closing of lines in NN-scattering and to the nucleon propagators within 
the AW reducible loop. It was found in [3] that those contributions cancel mutually. The 
basic mechanism is depicted in fig. 7 in terms of a one-loop approximation. By carefully 



Figure 7. One loop approximation in diagrams 5 (left) and 4 (right). The pion scatters in- 
side/outside the two-nucleon loop. 

considering the sums of all isospin and spin states one finds that both amplitudes just 
differ in the position of the derivative due to the squared nucleon propagators; for explicit 
expressions see [3]. After performing an integration by parts an exact cancellation of 
both contributions results. The process of mutual cancellation can easily be generalized 
to any number of two-nucleon reducible loops in diagrams 4 and 5 in fig. 6. Ann+1 



Figure 8. The cross indicates the action of the derivative with respect to k® . In this way, the first 
diagram on the second row of the figure is the same as the one to the left of the first row but with 
opposite sign, such that they cancel each other. The same applies to the second diagrams on the 
first and second rows. 

iterated wiggly line exchange in these figures implies n two-nucleon reducible loops. 
The scalar source can be attached for S5 to any of them, while for S4 the derivative with 
respect to k\, after performing the integration by parts, also acts on any of the iterated 
loops. Hence, 





s 4 + s 5 = . 



(13) 



Note, that this is a general argument following from the power counting and is not a 
property of the specific non-perturbative methods employed. 



Putting together the results of all contributions arising from the topologies in fig. 6 the 
ratio between the in-medium and vacuum quark condensate to NLO reads 
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(14) 



where the second and third term on the right hand side of eq.(14) rise from diagram 1 and 
the last term, which represents the quark mass derivative of the nuclear matter energy 
density, is identified with diagrams 3 and 6. This result agrees with the Hellmann- 
Feynman theorem for nuclear matter. Our findings, summarized in eq.(14), actually 
explain observations made in earlier calculations [16, 17] about the dominant role of 
the long-range part of nucleon interactions for the in-medium chiral quark condensate. 
Due to the cancellation of diagrams 4 and 5 and the pion-mass squared derivative of Sj, 
the pure the short-range part contributions vanish. 

IN-MEDIUM PION SELF-ENERGY 
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Figure 9. Contributions IT, to the in-medium pion self-energy up to &{p 5 ) entering the pion 
propagator A Jt (q,E > ) = 1/ (q 2 — + U(q,%) + is) in the nuclear medium. The squares corre- 
spond to NLO pion-nucleon vertices. For diagrams 3a, 3b, 8 and 10 the pion lines can leave or 
enter the diagrams. 

Here, we apply the chiral counting given in eq.(2) to calculate the pion self-energy in 
the nuclear medium up to NLO or 0(p 5 ), with the different contributions shown in fig. 9. 
We only consider here those NLO contributions to the pion self-energy in the nuclear 
medium that involve the NN interactions, for the rest see [1, 2]. They are depicted in 
the diagrams of the last two rows in fig. 9, where the ellipsis indicate the iteration of 
the two-nucleon reducible loops. Since V p = 2 in these contributions one needs the NN 
scattering amplitude at &(p°) to match with our required precision at NLO. In the same 
way as described for the case of the in-medium chiral quark condensate, we find a mutual 
cancellation of the diagrams 7 and 9 and of the isovector parts of diagrams 8 and 10 



[1,18] 



n 7 +n 9 = o, and n' 8 v +n\ v = o (15) 

Note that rig + riio also contains isoscalar 1/m recoil corrections, which are higher 
order, that do not cancel and yield contributions from NN interactions. 

CONCLUSIONS 

We have reviewed the development in [1, 2, 3] of an effective field theory for nuclear 
matter that consistently treats local short-range and pion-mediated long-range multi- 
nucleon interactions consisting of a novel in-medium power counting and the non- 
perturbative methods known from U^PT. We have calculated the nuclear matter energy 
per particle and found a naturally emerging saturation mechanism for symmetric nuclear 
matter and a weak parameter-dependence for neutron matter. For the in-medium chiral 
quark condensate we found that our results obtained from calculations of the generating 
functional agrees with the Hellmann-Feynman theorem. We also found the mechanism 
that explains why the short-range part of the NN-interaction contributes barely to the 
chiral quark condensate. For the in-medium pion self-energy up to NLO we find that all 
contributions from NN-interactions mutually cancel, which explains why pion-nucleon 
dynamics is phenomenologically successful for fitting pionic atom data. Meanwhile it 
has also been shown in [3] that analogous cancellations happen for the calculation of 
the in-medium pion decay at NLO. Those findings show, that there are no contributions 
violating the Gell-Mann-Oakes-Renner relation at NLO in the nuclear medium. 
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